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Abstract 

We present temperature dependences of the upper critical magnetic field H c2 and 
the Ginzburg-Landau parameter k for a ternary boride superconductor L^PdsB ob- 
tained from magnetization measurements. A specially developed scaling approach 
was used for the data analysis. The resulting H c2 (T) curve turns out to be surpris- 
ingly close to predictions of the BCS theory. The magnetic field penetration depth 
A, evaluated in this work, is in excellent agreement with recent muon-spin-rotation 
experiments. We consider this agreement as an important proof of the validity of 
our approach. 
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Magnetization (M) measurements in the mixed state of type-II superconduc- 
tors are traditionally used for evaluation of the upper critical field H c2 . If 
magnetic fields H higher than H c2 are available, the upper critical field is 
usually detirmined as a value of H, at which diamagnetic moment of the sam- 
ple vanishes. In real experimental situations, however, sample non-unformities 
and pining effects distort magnetization curves, making the choice of H c2 am- 
biguous. In this work, we propose and use a scaling approach to the analysis of 
magnetization data. As will be shown, this approach allows to obtain reliable 
H c2 (T) results even for samples of moderate quality with considerable pinning 
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effects. As an example, we consider magnetization data collected on a recently 
discovered ternary boride superconductor Li 2 Pd 3 B. 



The idea of our analysis is simple. In magnetic fields H much higher than the 
lower critical field H cl , the equilibrium magnetic susceptibility in the mixed 
state of type-II superonductors may be written as 



where the function f depends only on the ratio H/H c2 (T), i.e., the temper- 
ature dependence of x is determined by temperature variations of H c2 and 
k. The relation expressed by Eq. (1) is rather general. Although it follows 
directly from the conventional Ginzburg-Landau theory of the mixed state, 
this equation is not limited to conventional superconductors only (see [1] for 
details). 

Multiplying both parts of Eq. (1) by H/H c2 and introducing the universal 
magnetization function m = H f (H/ H c2 )/ H c2 , we rewrite Eq. (1) as 



■m (H/H c2 ) depends neither on temperature nor on the absolute value of k. 
This function may depend, however, on the symmetry of the order parameter 
or the sample geometry, i.e., it is not necessary the same for different super- 
conductors. For conventional superconductors with zero demagnetizing factor 
n, mo(H/ H c2 ) can be obtained, for instance, from numerical calculations of 
Brandt [2]. Although Eq. (2) is formally applicable only if H ^> H c i, the ac- 
tual limitations are not that strict. As may be seen in Fig. 3 of Ref. [3], even 
for k as low as 5, Eq. (2) represents a good approximation down to magnetic 
fields H ~ 0.1H c2 « 2H cl . 

Eq. (2) in a slightly different form was already used for the analysis of mag- 
netization data collected on high-T c compounds [1,4]. Here we slightly modify 
that procedure in oder to adjust it for materials with sufficiently low values of 
H c2 . In this case, the normal state paramagnetism may straightforwardly be 
evaluated from magnetization measurements above H c2 and there is no need 
to introduce an additional adjustable parameter as it was done in [1,4]. 



1 EXPERIMENTAL DETAILS 

The details of the sample preparation procedure can be found elsewhere [5]. 
The deviation of Li concentration from the stoichiometry was less than 1% for 



x(H,T) = f (H/H c2 )/K 2 , 



(1) 



m (H/H c2 ) = M(H,T)k 2 (T)/H c2 (T). 



(2) 



2 



this particular sample. A small (31 mg) part of the original (300mg) sample 
was used for magnetization measurements. 1 The sample had a platelike shape 
and it was placed with its longest side oriented along an applied magnetic 
field. The demagnetizing factor of the sample n = 0.12 was estimated by its 
magnetic susceptibility in the Meissner state. A commercial (Quantum Design) 
SQUID magnetometer was used for measurements. 



2 EXPERIMENTAL RESULTS AND ANALYSIS 
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Fig. 1. The field cooled magnetization measured in H = 5 Oe. The value of T c = 6.87 
K is estimated by the extrapolation of the H C 2(T) curve to H C 2 = 0, as will be 
explained in the text. 

Fig. 1 shows the field-cooled M{T) curve measured in a low applied field of 5 
Oe. As may be seen, the superconducting transition is rather broad (~ 10% 
of T c ), implying that the sample is not particular uniform. We also note that 
4irM(T — > 0) ~ 0.1H <C H, i.e., pinning effects are strong. A similar curve 
for the original sample is also shown. 1 Comparing these two curves, we may 
conclude that the smaller sample has substantially narrower transition and 
considerably weaker pinning. We also note that the smaller sample, studied in 
this work, has somewhat lower critical temperature. 

Isothermal magnetization curves measured in increasing (zero-field-cooled) 
and decreasing fields at T = 4 K are shown in the mainframe of Fig. 2(a). 
The M(H) curve measured above T c is shown in the inset. As may be seen, a 
normal-state paramagnetic contribution to M is substantial and it is not ex- 
actly linear on H. Nevertheless, because magnetization measurements at each 

1 The original sample (300 mg) was used for muon-spin-rotation measurements in 
[6] 
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Fig. 2. (a) M{H) curves below (4K) and above T c are shown in the mainframe and 
the inset, respectively. Extrapolation of M n to H < H C 2 is shown by the solid line, 
(b) The diamagnetic response of the mixed state. The solid line is the best fit of 
theoretical calculations of Brandt [2] to experimental data. The data points, which 
were used for fitting, are shown by black symbols. 

temperature were extended to magnetic fields well above H C 2, the normal- 
state magnetization curves M n (H), measured above H c2 , could reliably be 
extrapolated to H < H c2 , as it is indicated by the solid line in Fig. 2(a). 

The difference M sc = M — M n (diamagnetic response of the mixed state) is 
plotted in Fig. 2(b). In magnetic fields H > 14 kOe, M SC (H) is reversible, i.e., 
represents the equilibrium magnetization of the sample. The rounding of the 
experimental M SC (H) curve in higher magnetic fields, which may clearly be 
seen in Fig. 2, is due to sample non- uniformity and is expected in this kind 
of samples. In spite of this rounding, there is a magnetic field range (14 < 
H < 17 kOe), in which M sc is a practically linear function of H, as it follows 
from the Abrikosov theory of the mixed state [7]. In the following we assume 
that this is the behavior of the main part of the sample. Using Eq. (2), the 
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theoretical mo(H/H c2 ) curve, calculated in [3], can be fitted to experimental 
data with two fit-parameters H c2 and re. 2 This kind of fitting was applied 
to magnetization curves measured in the temperature range 1.8 K< T < 5.6 
K. The resulting temperature dependences of H c2 and re are plotted in Fig. 3 
by open circles. Because of stronger rounding of the magnetization curves at 
higher temperatures, this kind of fitting could only be used at T < 5 K. At 
higher temperatures, as may be seen in Fig. 3, the values of re are unreliable. 

Both H c2 (T) and re(T) agree rather well with theoretical calculations by Helfand 
and Werthamer (HW) [8]. This is why in the following analysis we shall use 
re(T) as it is provided by this theory (the solid line in Fig. 3(b)). In this case, 
we have only one adjustable parameter H c2 in Eq. (2), which greatly facilitates 
the fitting procedure and allows to extend the H c2 (T) curve to higher temper- 
atures. In all cases, we used for the analysis the steepest reversible part of the 
M(H) curves, as it is illustrated in Figs. 2 (b) and 4(a). This kind of fit worked 
quite well from the lowest available temperature of 1.8 K up to T m 6.2 K. At 
higher temperatures, no reasonable fitting to experimental data points could 
be achieved (see Fig. 4(b)). This breakdown is the obvious consequence of the 
fact that T c is not uniform across the sample volume. The rounded part of the 
transition becomes broader with increasing temperature and at temperatures 
close to T c no linear part of the M(H) curve can be identified (see Fig. 4(b)). 

The resulting H c2 (T) dependence is shown in Fig. 3(a) by closed squares. It 
is in perfect agreement with the theory in the entire available temperature 
range. By fitting of the theoretical H c2 (T/T c ) curve to experimental data, we 
obtain T c = 6.87 K and H c2 (T = 0) = 34.8 kOe. The value of T c , estimated 
in such a way, is consistent with low-field magnetization measurements (see 
Fig. 1). As has already been pointed out, the H c2 (T) curve is evaluated quite 
reliably, while the values of re, presented in Fig. 3(b), are overestimated. 2 

In Ref. [9], H c2 was evaluated as a point of the onset of the diamagnetism in 
M(H) and M(T) measurements. This is the reason that the value of H c2 (0), 
evaluated in [9], is higher than our result. We argue that the results of [9] do 
not represent the upper critical field as it is defined in the Ginzburg-Landau 
theory. 

We also show in Fig. 3(a) the H c2 (T) curve for Li 2 Pd3B obtained from magne- 
tization measurements by Doria et al. [10]. As may be seen, these data can also 

2 In the evaluation of re, we assume that the total volume Vq of the sample con- 
tributes to dM/dH. In nonuniform samples this is not correct. The effective volume 
V e ff, which should be used in calculations, is smaller than Vq. This means that our 
estimation correspond to re e // = Ky/Vo/V e ff > re. We also note that one should 
expect a decrease of 14// with increasing temperature, which may lead to some dis- 
tortion of the re(T) curve. However, because re is proportional to y/V, these effects 
are not expected to be significant. 
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Fig. 3. H C 2(T) and k(T) evaluated by fitting of magnetization data to the theoretical 
mo(i? '/ H C 2) curve with two fit-parameters H C 2 and k. The solid lines are the curves 
corresponding the Helfand-Werthamer theory [8] 

be very well fitted with the HW theory [8]. The resulting value of H&fi) = 35.7 
kOe practically coincides with our result cited above. We note that in Ref. [10], 
a very similar approach to the analysis of magnetization data was used. The 
only difference that the authors postulated that k is temperature independent. 
We argue, however, that some decrease of k with increasing temperature may 
clearly be seen as an enhancement of slopes of scaled magnetization curves 
with increasing temperature (see Fig. 1 in [10]). 

Using H c2 {T) and k(T), as they are presented in Fig. 3(a), we can scale our 
M(H) data according to Eq. (2) in order to obtain the m (H / H c2 ) function. 
The results are presented in Fig. 5. In complete agreement with our assump- 
tions, the data-points immediately above the irreversibility field collapse onto 
the theoretical curve. We also point out that at T < 4 K all the data, including 
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the rounded parts of the M(H) curves, merge together. Noticeable deviations 
develop only at temperatures closer to T c . A very similar behavior was earlier 
observed in NbSe2 samples [11]. 

Finally, we present the M(T) curve measured in a fixed magnetic field H = 10 
kOe (Fig. 6). Because both if c2 (T) and k(T) are already known, the theoretical 
curve can be calculated without any adjustable parameter. As may be seen 
in Fig. 6, agreement is rather good especially if errors in the extrapolation 
of M n (T) are taken into account. Here, at T > 5.2 K, we can also see some 
rounding of the experimental M(T) curve similar to that for isothermal M(H) 
data (Figs. 2 and 4). Considering the results presented in Fig. 6, we want 
to underline that the theoretical M(T) dependence is linear only in close 
vicinity of T C (H). This is why the linear extrapolation of experimental M(T) 
curves usually results in overestimated values of H c2 (see Fig. 6). Although this 
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Fig. 4. (a) Magnetization curves, (a) T = 1.8 K. The solid line represents theo- 
retical M(H) dependence fitted to experimental data as explained in the text, (b) 
T = 6.5 K. As may be seen, the theoretical curve (the solid line) cannot be fitted 
to experimental data. 
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Fig. 5. too = M(H,T)k 2 (T)/H C 2 versus H/H C 2- Only reversible magnetization data 
are shown. The solid line represents the theoretical dependence according to [2]. 

overestimation is not particularly strong, it may distort the resulting H c2 (T) 
curve providing a positive curvature of H c2 (T) at temperatures approaching 
T 

± c- 

We also note a rather unusual onset of irreversibility, which may be seen in 
the isothermal M(H) curves measured at temperatures well below T c (Figs. 
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Fig. 6. The magnetization M(T) curve measured in a fixed magnetic field of 10 
kOe. The solid line in the mainframe correspond to theoretical calculations for the 
Abrikosov vortex lattice [2]. The inset shows original M(T) data. The solid line 
represents the normal-state magnetization M n (T) extrapolated to T < T c assuming 
M n (T) ~ 1/T. 
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2(b) and 4(a)). Indeed, both brunches (increasing and decreasing fields) of 
the curves deviate to the left from the corresponding solid lines. This effect 
is non-observable at temperatures close to T c (Fig. 4(b)). The origin of such 
a behavior is unclear. Similar observations in high-T c superconductors are 
usually explained by a vortex-lattice-melting transition. However, it is difficult 
to imagine that this kind of transitions can exist in samples of such poor 
quality. 



3 DISCUSSION 



In the analysis above, we assumed that the temperature dependence of k 
follows the HW theory [8]. In complete agreement with this assumption, the 
slopes of the experimental M(H) curves vary with temperature as l//t 2 (T) in 
a rather extended range of temperatures 1.8 < T < 6.2 K (see Figs. 2(b), 4(a) 
and 5). We consider this fact as evidence of the validity of our scaling approach. 
At temperatures T > 6.2 K, sample imperfections distort the magnetization 
curves in a way that they cannot and should not be used for evaluation of 
H c2 (T). 

In the HW theory [8], H c2 (T) = ^2k(T)H c (T) (H c is the thermodynamic 
critical field) and both H C (T) and k(T) are calculated in the framework of the 
weak coupling BCS theory. In principle, our result, presented in Fig. 5, may be 
considered as an indication that the zero-temperature energy gap A(0) is close 
to the corresponding weak coupling value A(0) = 3.52T C . However, in order to 
be sure that this is indeed the case, both H c2 (T) and k(T) must coincide with 
theoretical predictions. In our case, as may be seen in Fig. 3(b), considerable 
pinning effects limit the accuracy to a degree that no unambiguous conclusion 
about temperature dependence of k can be made. In this sense, our results 
cannot be considered as contradicting to the value A(0) 4T c obtained for 
this compound in recent specific heat and muon-rotation experiments (/xSR) 
[6,12]. 

Now we compare our results with those obtained in /xSR experiments [6]. We 
remind that in this study we used a small part of the sample studied in [6]. 1 
The /xSR experiments provide the variance of the magnetic induction in the 
sample, which is inversely proportional to A 4 . Fitting of of these results to the 
Ginzburg-Landau theory gives H%\0) = 36.6 kOe and A^(0) = 2520 A[6]. 
We use an index (//) to denote the results obtained in /xSR experiments. 

Our value of H c2 (0) = 34.8 kOe practically coincides with H^(0). Taking 
into account that the two values are obtained from completely different ex- 
perimental data, this agreement should be considered as a rather convincing 
evidence that the Ginzburg-Landau theory may indeed be used for a quanti- 
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tative description of the mixed state even at temperatures well below T c if the 
temperature dependence of k is taken into account. 

Using H c2 (T), we may calculate £(T) = ^ /27rH c2 (T) ($ is the magnetic 
flux quantum). Extrapolation to T = gives £(0) = 97 A. Because k(T) is 
also known, we can evaluate the magnetic field penetration depth as A(T) = 
k(T)£(T). This gives A(0) = 3000 A. This value is only 20 % larger than 
A^(0). This is rather good agreement especially if we take into account that 
the samples are different. Furthermore, the two values of A(0) are not just 
close but the difference between them has a sign, which is expected. As was 
pointed out earlier, 2 our analysis gives slightly overestimated values of k and 
A. At the same time, as it is discussed in [6], in non-unuform samples, yuSR 
experiments provide A, which is smaller than its averaged value. 

In conclusion, using a recently discovered Li 2 Pd 3 B superconductor, we demon- 
strate that the scaling analysis of the equilibrium magnetization may success- 
fully be used in order to evaluate H c2 (T), k(T), and T c from reversible mag- 
netization data. It turned out that H c2 (T) is very close to predictions of the 
HW theory clearly indicating that Li 2 Pd 3 B should be considered as a con- 
ventional single-gap superconductor. Our results agree well with recent muon 
experiments providing H c2 (0) = 34.8 kOe, «(0) = 31.2, £(0) = 97 A, and 
A(0) = 3000 A. 
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